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^ We study the efTect of a finite quarlc mass on the quarlc number susceptibility in the 

framework of holographic QCD. We work in a bottom-up model with a deformed AdS black 
hole and D3/D7 model to calculate the quark number susceptibility at finite temperature 
with/without a finite quark chemical potential. As expected the finite quark mass suppresses 
' the quark number susceptibility. We find that at high temperatures T > 600 MeV the quark 

number susceptibility of light quarks and heavy quarks are almost equal in the bottom-up 
model. This indicates that the heavy quark like charm contribution to thermodynamics of a 
' QCD-like system may start to become significant at temperatures T ~ 600 MeV. In D3/D7 

model, we focus on the competition between the quark chemical potential, which enhances 
the quark number susceptibility, and the quark mass that suppresses the susceptibility. We 
'■ observe that depending on the relative values of the quark mass and the quark chemical 

(-H ' potential, the quark number susceptibility shows a diverging or converging behavior. We 

also calculate the chiral susceptibility in D3/D7 model to support the observation made 
with the quark number susceptibility. 

(N 

. §1. Introduction 

■ Fluctuations of conserved charges such as baryon number, electric charge, and 

strangeness are generally considered as useful probes for the structure of a thermal 
I medium, quark gluon plasma (QGP), produced in heavy ion collisions. Quark num- 

ber fluctuations are basic observables which can be obtained by taking derivatives 
of the grand canonical potential with respect to the quark chemical potential. The 
existence of a peak in the second order quark number fluctuation, or quark number 
susceptibility, near Tc is confirmed by recent lattice QCD calculations based on the 
Taylor expansion with respect to the quark (or baryon) chemical potential.- This 
' implies the existence of the critical end point (CEP), at which the first order phase 

transition terminates in the {fig, T) plane of the QCD phase diagram (see^* for a 
review). There have been many studies to calculate the quark number susceptibility 
and higher order quark number fluctuations in various model studie^"^ and lattice 
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simulations.'^'^ In|I3"ti21i quark number susceptibility at finite temperature 
with finite chemical potential is studied in a holographic QCD model. 

The AdS/CFT correspondenc^^SHllJi jg ^ powerful tool to study strongly coupled 
gauge theories. Using this correspondence, we can obtain physical quantities in gauge 
theories from gravity side. We can also discuss finite temperature systems from this 
gauge/gravity correspondence. In the hydrodynamic limit i.e. low frequency limit, 
we can calculate correlation functions.'^ The most fruitful hydrodynamic result 
is the observation of the "universal" bound of the ratio of the shear viscosity to 
entropy density.'^ Although the correspondence between QCD and gravity theory 
is not known, we can obtain much insights by using this correspondence. 

In this paper we study the effect of a finite quark mass on the quark number 
susceptibility in holographic QCD. We consider light, strange and heavy (charm) 
quarks. In the previous studies,'^''^ the mass of quark is zero, and so the quark 
number susceptibility with two massless flavor was calculated in a QCD-like system. 
We now extend the previous work by considering no-zero quark masses. Our basic 
motivation for the charm quark is an observation made in perturbative QCD anal- 
ysis on the role of quark mass in QCD thermodynamics. In^^ it is claimed that 
unexpectedly, the charm quark plays some role in the QGP at relatively low temper- 
atures T ~ 350 MeV. Since this analysis is based on the perturbation, it should be 
important to check the claim in a non-perturbative study. To explore the effect of 
a finite quark mass at finite temperature and/or density, we work with two models. 
The first one is a bottom-up type model based on the deformed AdS black hole due 
to finite quark mass obtained inW^ Then we consider D3/D7 model with a finite 
chemical potential studied in.'23'ESl' 

This paper is organized as follows: After introducing the quark number suscep- 
tibility, we calculate it in the deformed AdS black hole background in Section 2. The 
hydrodynamic and thermodynamic analysis are given in the section. In Section 3, 
we discuss the quark number susceptibility in D3/D7 system. We summarize the 
result in the final section. The procedure to obtain the retarded Green function is 
briefly given in the Appendix A. In Appendix |Bl we consider the chiral susceptibility 
in D3/D7 model, though the chiral symmetry in D3/D7 is only U{1)a not that of 
QCD, SU{Nf)i^ X SU{Nf)n. 

§2. Quark number susceptibility in deformed AdS black hole 



The quark number susceptibility was proposed as a probe of the QCD chiral 
phase transition at zero chemical potential,'^''^' 

--If- 

Later it has been shown that the quark number susceptibility can be rewritten in 
terms of the retarded Green function through the fluctuation-dissipation theorem,'^ 



Xq = -limRe Gtt(a; = 0,/c) , (2-2) 
where y{oj,k) is the retarded Green function of the vector- vector correlations. 
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Now, we calculate the quark number susceptibility in the deformed AdS black 
hole background. We do this in two different ways that should be equivalent each 
other. We first use the equation (j2-2|) to conduct hydrodynamic analysis and then 
calculate the quark number susceptibility using the thermodynamic relation in the 
equation (j2-ip . 

The background we shall discuss is AdS black hole with back-reactions coming 
from flavors 

{dsf = e-2^(^) ( - f{z){dtf + {dxf) + (2-3) 



with 




The parameter zq is defined as 



(2-4a) 
(2.4b) 



Mq represents the quark mass and = SvrGs with 1/G^ = (i2N^ /{ttL^) where L is 
the AdS radius. IrP3 there were two deformations to be discussed, deformed AdS 
and deformed AdS black hole. In the present work we only consider one of them, 
i.e. the deformed AdS black hole background. Each of these two backgrounds is 
reasonable with given potential^. Note that irP^ Nf counts the number of massive 
quarks only. For instance, if we are interested in a 2+1 flavor system, we take Nf = 1 
in the equation (|2-5|) , neglecting the back-reaction from the two flavors whose masses 
are zero or very small. Taking the limit Mg — t- 0, the metric (|2.3p is reduced to the 
AdS black hole 

(d.)2 = l(-/(z)(dt)2 + (dx)2)+ 



z^izY 



with 

4 



= 1 



where Zh implies the location of horizon of the original AdS black hole. There exists a 



horizon zh which gives zero of the function ()2.4b[i**)| . Then the Hawking temperature 



*) ijPZ) tiig authors discussed the Hawking-Page type transitiorP^''^ which is the dual to the 
confinement/deconfinement transition using these two deformed backgrounds. However, since the 
bulk potentials for each background are not the same, the Hawking-Page type analysis needs 
to be improved. 

**' We may refer the position of horizon as 



Zh ~ ZQ \ i + ProductLog 



{zh/zQ)*-2 



2e 

where ProductLog(a) stands for a solution for x in the equation xe^ — a. 
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e*^"«^ . (2-6) 



is given by 

■nzH V Zh 

We note that for small Mg, zh becomes Zh- Using the background ()2-3p . we assume 
Mq is not large within the small mass deformation. Since the background (|2-3p is 
constructed from the bottom-up type approach, we could not find any restrictions 
for the finite value of Mq. 

Let us now consider 5D U{\) gauge field which is dual to 4D quark number 
current on the background ()2-3p . We then calculate the quark number susceptibility 
by using the AdS/CFT correspondence. Since in this background there is no back- 
ground charge which could be interpreted as the chemical potential in the boundary 
theory, we here discuss the mass dependence of the quark number susceptibility. 

We start by introducing a dimensionless coordinate u = {z/zhY which is nor- 
malized by the horizon. In this coordinate system, the horizon and the boundary 
are located at n = 1 and u = 0, respectively. The metric (j2-3p is rewritten as 

{dsf = e-2^(") { - f{u){Atf + {dxf) + (2-7) 



with 



K{u) = \\og{zlu) + \[^)\, (2.8a) 




The action for the U{1) gauge field is 

1 



(2-8b) 



S = --^ / d^x^FmnF'^'', (2-9) 

where is the 5D gauge coupling constant. We shall work in Au{x) = gauge and 
use the Fourier decomposition 

A^{t,z,u) = I -0^ e-*^*+*'=M^(a;,A;,n), (2-10) 

where we choose the spacial momenta which are along the z-direction. Variations of 
the action with respect to At{u) and Au{u) give equations of motion 



. 1 
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0=(e A',{u)J -—jL_(^k^At{u)+u;kA,{u)y (2-lla) 

= ujA[{u) + kf{u)A'^{u), (2-llb) 
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where the prime stands for the derivative with respect to u. The equation of motion 
for Az{u) can be derived from the equations (j2-llap and ()2-llbp . For Ax{u) and 
Ay{u), one can obtain decoupled equations of motion. Since we are interested in the 
time-time component of the retarded Green function to calculate the quark number 
susceptibility (|2-2p . we will not consider Ax{u) and Ay{u) hereafter. 

From the equations (|2-llap and (|2-llbp . we obtain an equation for At{u), 



[uf{u) e"('4)"4'(n))' 



2\zqJ y ' J 4 \f{u) 



uf{u) e--^'^> ") - ^ e"^^^^ "( ^ - A:^) \A[{u). 

(2-12) 



2 



We can make the equation more convenient form by changing the variable, 

A't{u)=e^^'Q'''x{u). (2-13) 
Then the equation ()2T2p can be rewritten as 

= {uf{u)X'{u))' + «'^^^'"(;^ - ^')^(^)- (2-14) 

Let us proceed to solve the equation of motion ()2-14p . This is an ordinary second 
order differential equation with a regular singularity at the horizon u = 1. Writing 
a solution as X{u) = (1 — uYF{u) where F{u) is a regular function at the horizon, 
we fix the constant v by imposing the incoming wave condition, 

T^ = -i-r^, 2-15 



where T is the Hawking temperature ()2-6p . 

Now we shall solve the equation of motion in the hydrodynamic regime i.e. small 
uj and k compared with the temperature T. In this regime, we could expand the 
function F{u) as 

F{u) = Fq{u) + ujF^{u) + k^ Fk2{u) + 0{uj'^ ,u:k^). (2-16) 

Since we are interested in the quark number susceptibility defined by ()2-2p . we only 
need solutions for Fq{u) and Fj^2{u). Relevant equations are read off as 

0= (n/(n)F^(n))', (2-17a) 



= [uf{u)F',,{u)) - '-f e^^^Q^ >oM. (2-17b) 

In the equation ()2-17ap . avoiding the singularity at the horizon, we fix the function 
Fq{u) as a constant C which is determined later, 

Fo(n) = C. (2-18) 



2 1 ('H 
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Using this solution, the equation (|2-17bp leads 

Auf{u)J^ Auf{u) \zhJ \ y' 

(2-19) 

where we have imposed the regularity at the horizon. We insert the solutions above 
into the equation ()2-llap and obtain 

2 

-^^IL_(^k^At{u) + LokA,{u)^ =X'{u) = eFl,{u) + 0{u:). (2-20) 

Taking the definition of the quark number susceptibility (j2-2p into account, the 
equation (|2-20p is enough to determine At{u), 

= ^L_^ (f _ , (2-21) 

where we have defined the boundary variable as = At{u)\u~^Q to fix the constant 
C. 

Let us evaluate the quark number susceptibility. The on-shell action to be needed 
is estimated as 



u=0 

(2-22) 

u=l 



Following the procedure given in Appendix A, we can read off the two point function 
Gt t(u;, k). We then conclude that the quark number susceptibility is given as 



55 

In Fig{Tl we plot Xg/^^ ^ function of T with varying the quark mass mJ*^ . 
Note here that the quark mass in our model can be different from that in QCD (or 
in the chiral Lagrangian) by a constant factor. For instance, in the hard wall model 
the factor is \/3i which is obtained by matching the scalar two-point correlators 
calculated from QCD and the hard wall modelP M^"^^"^ = Mq/2,. If we assume 
that this is valid in our case, we have to divide the quark mass in our model by the 
factor \/3 when we compare our results with those from QCD or chiral perturbation 
theory (ChPT). In Fig{T| the vertical dashed line is for a typical deconfinement 
temperature with no quark mass effect in the hard wall model.'^S If we assume 
that the deconfinement temperature in our approach is roughly equal to that in 



*-* We are interested in the qualitative tendency of the mass dependence to the quark number 
susceptibility. For the illustration purpose, we vary Mq from 100 to 2000 MeV. 
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Fig. 1. The quark number susceptibility as a function of the temperature for various quark masses. 
The vertical dashed line denotes the deconfinement temperature in the hard wall model. 



the hard wah modeO, our resuhs on the left of this hne are just for illustration 
and might be replaced by those with the thermal AdS background as in.lI^''ISl' As 
expected the finite quark mass suppresses the quark number susceptibility. We find 
that at high temperatures T > 600 MeV the quark number susceptibility of light 
quarks and heavy quarks are almost equal. This indicates that the heavy quark like 
charm contribution to thermodynamics of a QCD-like system may start to become 
significant at temperatures T ~ 600 MeV. Though it might not mean much, we 
compare our number T ~ 600 with that from lattice QCD. Inj^S' it is shown that 
the fluctuations of the charm quark number is comparable to those of the light 
quark sector at high temperatures and claim that at temperatures T ~ 800 MeV the 
contribution of charm quark to QCD thermodynamics is significant. 

Before closing this section, we briefiy demonstrate the thermodynamic analysis 
on this system. In the background ()2-7p . we turn on the temporal component of 
gauge field which could introduce the finite density in the boundary gauge theory. 
We work in the grand canonical ensemble. In the AdS/CFT correspondence, we 
could identify the grand potential i?(T, /i) as 

17 = T5on-shell) (2-24) 

where S'on-sheii is the 5D bulk on-shell action evaluated in Euclidean spacetime. As in 
the usual imaginary time analysis, we analytically continue to Euclidean spacetime 
and compactify the time direction with the period 1/T. Assuming that the gauge 
field depends only on the radial coordinate u, the equation of motion leads 

A;(n) = -d e^^^^ (2-25) 

where the integration constant d (> 0) will be proportional to the quark number 

*' We should consider the Hawking-Page type transitiorP^''^ in the current setup. Since the 
analysis for this confinement/deconfinement irP^ needs to be improved as we explained, we here 
take the hard wall model as the closest existing one. 
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density. Then the grand potential ()2-24p can be estimated as 

72t/ ^2 



n = ^ due'^^Q' = /. '^ (l-e^^-Q^ ), (2-26) 



where V3 is 3D vohime. The chemical potential is defined by 



2dzf, , - , — , > 
^1= I duA't(u) = j^{l-e^^'Q' ). (2-27) 



Using these relations, we obtain 



Since the quark number density is defined by 

1 df2 



■d), (2-29) 



one can easily confirm that the grand potential exactly gives the quark number sus- 
ceptibility obtained in ()2-23p . As it should be, the hydrodynamic and thermodynamic 
analysis give the same result. 



§3. Quark number susceptibility in D3/D7 system 



In this section, we study effects of the chemical potential and the quark mass 
for the quark number susceptibility at finite temperature in a D3/D7 system. Fol- 
lowing we consider the "black hole embedding" of the probe D7-branes in the 
D3-brane background in which we could introduce the finite quark mass together 
with the chemical potential. In this section we follow the procedure developed iiJ23 
and then extract relevant parts to discuss the quark number susceptibility. 

We start by introducing Nc D3-branes. The AdS/CFT correspondence provides 
the dual descriptions of = 4 SU{Nc) SYM as the type IIB string theory on 
AdSsxS^ with the identifications L^/lg = '^gyu^c = 2A and Qs = (7ym/(2^)- Here 
L, Is = ^/oH and A are the radius in AdSs and S^, the string length and 't Hooft 
coupling, respectively. Taking the large and the large 't Hooft coupling limit, 
the string description reduces to that for the supergravity. The finite temperature 
system in the deconfined phase can be realized by using a corresponding black hole 
geometry in the Euclidean signature: 

(d^)' = ^ [hm^ + (d-l^) + 5 + v'^^nt) , (3-1) 

with 

/o(-) = l-l 
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where uq implies the location of horizon. The periodic thermal identification t 
t + 1/T leads the Hawking temperature 



T = ^° 
7rL2- 

Following the paper p3 we introduce new coordinate q through 



(3-2) 





and rewrite the metric (13-11) as 



with 



= u^ + Ju'' -ul, (3-3) 



(d-)' = \il i £(dt)^ + fm') + G'dnl), (3-4) 



We work with the dimensionless coordinate p = q/uq in the lest of the paper. We 
refer to the horizon as p = 1 and the AdS boundary as p (= V^u/uq) — oo. 

In order to introduce the fundamental matters, we may consider flavor branesp^ 
which provide flavor gauge fields in the world volume. Here we put Nj D7-branes 
with the following intersection 

0123456789 
D3: X X X X (3-5) 
D7: xxxxxxxx 

where D7-branes are wrapping on S'^ of S^. The additional degrees of freedom which 
we would like to discuss are generated by open string oscillations between D3 and 
D7-branes described by A/" = 2 hypermultiplet. The fermions in this multiplet could 
be interpreted as the analog of quarks in QCD. Taking the probe approximation 
Nf <^ Nc, we could consider the overall dynamics of Nf D7-branes which can be 
described by the DBI action in the lOD background 



S = NfTr j dV e-^ V det(5,nn + 2^a'Fmn), (3-6) 

where gmnicr) and Fmn{(y) are the induced metric and the field strength of the world 
volume gauge field, respectively. The tension of D7-branes and the dilaton field 
are given by Tj = l/((27r)'^/f ) and e''^ = Qs, respectively. It is convenient to divide 
the transverse 6D part to D3-branes in ()3-4p into two parts i.e. 4D and 2D whose 
coordinates are given by spherical (r, Q^) and polar (ii, (p) coordinates, respectively, 

{dpf + p'^dQl = {drf + r'^dQl + {dRf + R^{d^f 

= {dpf + p2 ({def + sin^ 9 dni + cos^ 6 {d(p)^) , (3-7) 



where r = psinO, R = pcosO with < < 7r/2 and p^ = r"^ + R^. By construction, 
we consider the case where the world volume coordinates of D7-branes are given in 
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the static gauge i.e. o"™" = {t,x,p, H^). Due to the symmetries for the translation in 
{t, x) and the rotation in (p, ]?3), the embedding of the D7-branes could only depend 
on the radial coordinate p. Since the rotational symmetry in {R, ip) allows to be 
(p = 0, the embedding might be characterized by xip) = cos^ through 6{p) which 
is the angle between two spaces (r, /^s) and {R,if). The asymptotic value of the 
distance between D3 and D7-branes which is given by R{p) for large p provides the 
quark mass Mg. The induced metric on the D7-branes is given by 



(ds)D7 = L I -^p y-jidt) + /(dx) j + -2 

+ (1-X') dr?||, (3-8) 

where the dot stands for the derivative with respective to p. We also introduce the 
non-dynamical temporal component of the gauge field At{p) which leads the chemical 
potential and the density at the AdS boundary. 

By using the induced metric ()3-8p and the gauge potential At{p), the DBI action 
(13-61) now becomes 



So = ^^'^f j dp - - X^+P^X^) - 2/(1 - X^)Al (3-9) 

where we have defined At{p) = 2t^ a' At{p) /uq. Since there exist no At{p) terms in 
the action, the equation of motion for At{p) can be reduced to the following form 
with an integration constant (i, 

J. '^P^-^'i'^' (340) 



2V/'(1-x' + pI\")-2/(1-xM 

The equation of motion for x{p) is given as 



1/2- 



P^ffX j /o.-, _2^ , . 2.2\/-, , 



1/2 



+ . "V" ... 3(1 - x') + 2p^x' 1 + 



-1/2- 



where we have eliminated the gauge field At{p) by using the relation (|3-10|) . Near the 
boundary, asymptotic solutions of the equations of motion ()3-10p and ()3-lip behave 
as 

It(p) = /I-4 + --- , (3-12a) 
P^ 
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x(/>) = - + 4 + --- • (3-12b) 

P P-^ 

The chemical potential fi can be defined as the boundary value of At{p), while the 
quark mass Mg can be estimated through the asymptotic value of the separation 
of D3 and D7-branes i.e. Mg = px{p)/{2iTa') at the AdS boundary. Taking the 
rescaling of the gauge field At and the coordinate p into account, the integration 
constants p and m can be related with these values, 

~ 2W [2p 

p = p = \ --, 3-13a 

no VAT 

m = 27Ta'—M„ = ^^. (3-13b) 

As we will estimate belo'wF^. the remaining constants d and c would be proportional 
to the the quark number density Uq and the quark condensate (tpij), respectively. 

We have to solve the equations of motion ()3-10p and ()3Tip numerically to de- 
termine the embedding and the property of the gauge field. Here we restrict to the 
black hole embedding in which the D7-branes touch the horizon, since this might 
be thermodynamically favored configuration in the system with finite density. We 
impose boundary conditions at the horizon as x{l) = and At{l) = to remove 
singularities and = Xo- We fix m and p which depend on xo and d by matching 
the numerical solutions with the asymptotic forms at the boundary. 

We now consider the on-shell action which is related with the partition function 
Z of the field theory in the context of AdS/CFT correspondence: 

Z = e-^°"-*on. (3.14) 

However the on-shell action in this system contains UV divergences. It is well- 
known that one can prepare local boundary counter terms for probe D-branes in AdS 
spacetime by applying the holographic renormalization.'^ In the D3/D7 system, 
taking the asymptotic solution (|3-12b|l into account, the relevant boundary counter 
terms take the form,'321' 

Set = ^^g^^^V sj - ^((pLx - m'f - 4mc) |, (3-15) 

where Pmax is the cut-off for UV divergences which may go to infinity after precise 
calculations. It should be noticed that there exist finite contributions in the counter 
terms. Together with these counter terms, we could obtain the regularized action 

S = So + Set 
_ XNcNfT^ ( 



32 




M p'fil - x')^P{l -X^ + P^X^) - 2/(1 - X^)Aj 

-p^ + m'^f^ -^{{m^ -if - Ara(^\. (3-16) 



In the dimension three operator, we neglect contributions from squarks in the hypermultiplet. 
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Since we are interested in the black hole embedding, the integration supports from 
the horizon to the AdS boundary. Evaluating the on-shell action for (j3-16p which is 
reduced to boundary values through the equation of motion, we could observe that 
the quark condensate is proportional to the integration constant c, 

^^^^ = "^9^^°^^ = -^VAiVeiV/r^c, (3-17) 

where we have used the asymptotic solutions ()3-12ap and ()3-12bp and the relation 
()3-13bp . As we did in the previous section, we could identify the grand potential 
i7(^, r) as 

f2 = -TlogZ = TS'on-shell- (3-18) 

The quark number density can be also calculated through the on-shell evaluation, 

"« = -^|r4/|w^'^. (3-19) 

Let us now calculate the quark number susceptibility numerically using ()2-ip and 
()3-19p . In Figill we plot our results with a = Jl/m ~ fJ-fMq. Since 1/m ~ T/Mg, we 
can consider the horizontal axis of the figure as the temperature with the choice of Mg 
fixed. For very small a 1) or Mg S> fJ-fV^, the quark mass is dominating the quark 
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Fig. 2. The quark number susceptibility in the D3/D7 model as a function of 1/m. Here a 
Ji/m = fi/{V2Mq). 



number susceptibility and suppresses it at low temperature, while in the opposite 
case the chemical potential enhances it. Note that the confinement /deconfinement 
transition temperature in the D3/D7 model is zero, since there is no dimensionful 
parameter which could characterize the critical temperature as in the Hawking-Page 
transition.l^''^ Indeed, within the large Nc approximation, one can observe the 
entropy is zero for T = 0, while for given finite T that is proportional to N^T^ whose 
dependence indicates the deconfinement of the color degrees of freedom.l^D In 
addition, if one would like to consider mesons which correspond to the fluctuations 
of D7 flavor branes, one could observe their discrete mass spectra at T = oW^ 
Another cautionary remark is that in the QCD phase diagram, there exists a CEP 
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at a single value of the chemical potential, while Fig. 2 may imply multiple points in 
the sense that the quark number susceptibility is diverging for many different values 
of the chemical potential with fixed quark mass. At this moment, we have no clear 
understanding how to remove multiple CEPs. We may hope that some large Nc 
corrections could resolve this problem. Without resolving this issue, we may not be 
able to address single CEP in the QCD phase diagram within D3/D7 model. 

Now we discuss implications of our findings in QCD at high temperature. We 
first spell out some limitations of our approach. As well known, a generic problem 
with this kind of holographic QCD is that the results from holographic QCD are 
mostly large Nc leading ones. Though there are some studies to include the cor- 
rections, most of them include some subset of the corrections, not all of them in 
a consistent way. Apart from this generic plague, the D3/D7 model we used here 
shows that deconfinement temperature Tc is zero. In addition, we assumed an exact 
isospin symmetry to use the abelian DBI action for the Nf probe branes, and so our 
approach could not study A^'j-dependent nature of QCD phase transition. Unlike 
QCD, the temperature T, the chemical potential /i, and the quark mass Mq are not 
all independent, and only two of them or ratios of them are independent, which is the 
relic of conformal nature of our background metric. With the above-mentioned in 
mind, we comment on QCD, more precisely a QCD-like system, at finite temperature 
based on Figi2l 

For small a 1) or Mq » jJ-lV^, the quark number susceptibility does not 
show any diverging behavior as the temperature approaches Tc, indicating that if 
the quark mass is big enough compared to the value of the quark chemical potential, 
the transition to deconfined phase is a smooth crossover rather than first or second 
order phase transition. In the opposite case a > 1, it shows a rapid jump-up, coming 
close to Tc from high temperature. This implies that as long as the quark chemical 
potential is a few times larger than the quark mass, the transition would be second 
order. In Appendix |Bl we calculate the chiral susceptibility in the D3/D7 model, 
see Figini We observe that it shows a similar behavior with the quark number 
susceptibility in Figj2j 

§4. Summary 

We studied the effect of a finite quark mass on the quark number susceptibility in 
the framework of holographic QCD. We first considered the bottom-up model with 
the deformed AdS black hole that includes the back-reaction of the quark mass. 
As expected the finite quark mass suppressed the quark number susceptibility. We 
observed that at high temperatures T > 600 MeV the quark number susceptibility of 
light quarks and heavy quarks are almost equal, indicating that the heavy quark like 
charm contribution to thermodynamics of a QCD-like system may start to become 
significant at such temperatures. 

We then moved to the D3/D7 model to calculate the quark number susceptibility 
at finite temperature with a finite quark chemical potential. We studied the com- 
petition between the quark chemical potential, which enhances the quark number 
susceptibility, and the quark mass that suppresses the susceptibility. We observed 
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that depending on the relative values of the quark mass and the quark chemical 
potential, the quark number susceptibility shows diverging or converging behavior. 
With a caution that our approach is still not that close to a realistic system like 
quark-gluon plasma, we discussed physical implication of our results in a QCD-like 
system at the end of section [3l 

We also calculated the chiral susceptibility in D3/D7 model in Appendix iBl to 
support the observation made with the quark number susceptibility. 
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Appendix A 

Minkowskian correlators in AdS/CFT correspondence 

In this appendix, we briefly summarize the prescription for the Minkowskian 
correlator in AdS/CFT correspondence. We here follow the prescription proposed 
in.f^ We work on the following 5D background, 

{dsf = gf.Ax^dx'' + Quuiduf. (A-1) 

where x^^ and u are the 4D and radial coordinates, respectively. We refer the bound- 
ary at u = and the horizon u = 1. Let us consider a solution of an equation of 
motion in this 5D background. Suppose a solution of an equation of motion is given 
by 

<t>{n,x) =/(04 e^'Vfc(n)/(fc), (A-2) 

where fk{u) is normalized such that /fc(0) = 1 at the boundary. After putting the 
equation of motion back into the action, the on shell action might be reduced to 
surface terms 



s[cp'] = i ^,<iP{-k)g{k,u)<i^\k) 



M=0 



u=l 



(A-3) 



Here, the function Q{k,u) can be written in terms of f±k{u) and duf±k{u). Ac- 
commodating Gubser-Klebanov-Polyakov/Witten relatiorPD'ES iq Minkowski space- 
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time, Son and Starinets proposed the formula to get the retarded Green functions, 



G{k)=2g{k,u) 



(A-4) 



u=0 



where the incoming boundary condition at the horizon is imposed. In this paper, 
we consider correlators of U{1) currents J^{x), where Jfi{x) is the vector current of 
quark field or quark number current. Now we define the precise form of the retarded 
Green functions which we discuss in this paper: 



u{k) = -i / d^x e-^^^e(t)([J^(x), J,(0)] 



(A-5) 



Appendix B 

Chiral susceptibility in D3/D7 system 



In this appendix, we consider the chiral susceptibility. Note that in D3/D7 
system, we have only U{1) axial symmetry. Chiral susceptibility Xc is one of the 
important observables in terms of chiral symmetry restoration: 



Xc 



dM„ 



(B-1) 



In the D3/D7 system, the quark condensate (i/^ip) is given by ()3-17p . Note that 
our chiral condensate (V'V') is defined to be positive since c is negative, while it is 
negative in ordinary QCD. Reading out the numerical value of c, we calculate the 
chiral susceptibility. 
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Fig. 3. The chiral susceptibility as function of 1/m. 
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